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Abstract: In holographic models of gauge theories with matter, there generically exists 
a first order phase transition in which mesons dissociate. We perform a careful analysis of 
the meson and quasi-particle spectra in the overheated resp. undercooled regimes close to 
the junction of the two phases. We show that all overheated finite meson masses eventually 
become infinitely degenerate, which implies that any connection to the quasi-particle spec- 
trum is more subtle than previously suggested. For the Sakai-Sugimoto model no smooth 
connection between these spectra is possible. 
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1. Introduction 

Holographic models for gauge theories with matter exhibit a rather generic three-phase 
structure. For low temperatures, the background is a regular geometry, representing the 
confined phase of the gauge theory. At sufficiently high temperature, the background 
undergoes a Hawking-Page transition to a black brane geometry. In the dual gauge theory, 
this transition corresponds to the deconfinement of gluons. Mesons still remain bound, 
as long as the brane representing the matter degrees of freedom does not intersect the 
background horizon. When this finally happens, a second transition occurs, to a phase 
in which mesons melt and the gauge theory is fully deconfined. We will call these three 
phases the low, intermediate and high temperature phase respectively. 1 

After the mesons have melted, one is left with a strongly-coupled quark-gluon fluid. 
This fluid exhibits a discrete spectrum of unstable excitations [4], which are described by 
quasi-normal modes on the probe brane. It has been suggested [5, 6] that the thermal pole 

1 There is by now a large list of references dealing with this three-phase structure; seminal papers are [1] 
for the D3/D7 system, [2] for the D4/D6 case and [3] for the D4/D8 model. In special situations the 
intermediate phase can be absent, but we will not discuss those cases here. 
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masses of mesons before the transition should be related to the real parts of the quasi- 
normal modes after the transition. It would then be possible to follow the stable mesons of 
the intermediate temperature phase into the high-temperature phase, where they become 
unstable quasi-particles with finite decay width. The fact that the widths are of order one 
is consistent with the expected behaviour of meson decay widths in the deconfined phase. 2 
Moreover, lattice results indicate that finite- width mesonic resonances indeed survive above 
the deconfinement temperature (see e.g. [8, 9]). 

While such a connection between stable and unstable meson modes thus seems rea- 
sonable, it so far relies on numerical evidence. Other known facts about the generic ther- 
mal behaviour of meson masses suggest that the connection might be more subtle. In 
the intermediate temperature phase, the pole masses of mesons are known to decrease 
as the temperature goes up [10]. This shows qualitative resemblance with lattice QCD 
studies [11, 12]. In QCD proper, chiral partners such as the p and a\ meson eventually 
become degenerate when chiral symmetry is restored. In contrast, all holographic models 
go through the first-order meson melting phase transition before such degeneracies can be 
observed. However, the graphs of the thermal behaviour of the masses [10] suggest that 
a degeneracy may still occur when the system is overheated. If such a degeneracy occurs, 
then an identification of the discrete spectrum before and after the transition becomes 
more subtle (we will indeed see that this happens). 

A further problem occurs in the Sakai-Sugimoto model with zero bare quark mass. 
The embeddings after the transition are all equivalent, which implies that the quasi-normal 
mode spectrum after the transition is fixed, i.e. it never becomes degenerate. If the meson 
spectrum, on the other hand, does become degenerate at the junction, then no smooth 
identification seems possible. In addition, 3 mesons of this model transform under the 
adjoint of the diagonal SU(Nf), while the quasi-normal modes sit in the adjoint of SU(Nf)L 
or SU(Nf)n. This again makes a straightforward identification of the spectra unlikely. 

In the present paper we therefore try to answer two key questions: what is the precise 
behaviour of meson pole masses as the critical embedding is approached from below the 
transition temperature, and does this leave any room for a connection to the real parts of 
the frequencies of the quasi-normal modes just above the critical temperature? We will 
restrict ourselves to an analysis of the spectrum of vector excitations close to the critical 
embedding, and comment only briefly on a similar analysis for scalar excitations. For the 
sake of simplicity, we will also only consider modes homogeneous in the internal spheres; 
generalisation to other cases should be straightforward. 

We will first provide numerical evidence which indicates that the potential of the 
associated Schrodinger problem takes a particularly simple form as the critical embedding 
is approached. Inspired by these numerical results, we then present an analytic procedure 
to compute the spectrum close to either side of the phase transition. We will show that 

2 In the intermediate-temperature phase the decay of a meson to a quark-antiquark pair would also be of 
order one, but because mesons are anomalously light (of order 2m q /^/\ where m q is the constituent quark 
mass), there is no phase space available for their decay. Meson to two-meson decays [7] are of order 1/N C . 

3 We thank Ofer Aharony for emphasising this point to us. 
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any finite number of mesons in the discrete spectrum of the intermediate phase become 
degenerate in mass at the critical (overheated) embedding (our findings avoid the no- 
crossing theorem of Von Neumann & Wigner [13] because a finite change in the temperature 
leads to an infinite change of the potential). 

A similar degeneracy occurs, but only in the D3/D7 case and the polar embedding 
of the D4/D6 systems (see figure 0), for the quasi-normal frequencies when the critical 
embedding is approached from the high-temperature side (undercooked). Since the masses 
and quasi-normal modes are labelled by a single excitation number n, there is in principle an 
infinite number of ways to relate the states in these two spectra, and the connection between 
mesons and quasi-normal modes thus seems more subtle than previously suggested [5, 6]. 
For the D4/D8 and equatorial D4/D6 systems, the quasi-normal modes never degenerate, 
and there is thus no smooth connection to the meson spectrum. 



2. Setup and numerical motivation 

Before we discuss our main analytic results for the behaviour of the meson spectrum near 
the melting transition, let us first set the stage and motivate our approach by some nu- 
merical results. 

Our starting point is the generic expression for the background metric of a Dp-brane. 
It is given by 

d, 2 = (- f p {u)^ + <%d^) + {^y (y^ + u2d ^-p) • t 2 - 1 ) 

where i,j = l...p and f p (u) = 1 — (ut/u) 7 ~ p . There is also a non-trivial dilaton 

„v (p_7)(p-3)/4 



L 



) • (2-2) 



The temperature of this background is given by 

UT ^ (2.3) 



4vrL V L 

Probe branes can be embedded in various ways, depending on the value of p. In all cases 
which we will discuss, four of the directions of the probe brane are aligned with t and 
three of the x % . For the D4/D6 system one has a choice of two different subspaces for the 
embedding of the probe. One of these occurs again for the D3/D7 system and the other 
one occurs for the D4/D8 system. See figure [j] for details. 

There are two types of embeddings of the probe brane, depending on whether or not it 
intersects the horizon of the background. The "Minkowski embeddings" do not reach the 
horizon. The non-equatorial embeddings of this type are described most conveniently in 
the r-A coordinate system using the A coordinate on the world-volume. The "black hole 
embeddings", on the other hand, do reach the black hole horizon, intersecting it orthogo- 
nally. The non-equatorial embeddings of this type are best described in a polar coordinate 
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Figure 1: The generic holographic setup for the D4/D6 system. The r-A plane is a 



conformally flat plane arising from the second term in (2d), and r = x (see the appendix 



for details) . The embeddings in the horizontal r-A plane will be called "equatorial" . The 
embeddings of the D4/D8 system are similar to these equatorial ones. The D3/D7 system 
only has embeddings corresponding to the ones in the vertical ? — A plane, which we will 
call "polar" . Critical embeddings are displayed in red, while the horizon is indicated by 
the blue surface. 



system p-6, as indicated in figure [T] (see appendix A.l for details on our coordinate sys- 
tems). The embedding which lies along the horizontal axis plays a special role, as it occurs 
in all models which we analyse here. 

Statements concerning a possible relation between meson modes and quasi-normal 
modes have all been concerned with the behaviour of embeddings in the r-A plane, in other 
words, with the D3/D7 and D4/D6 system. There is, however, also a critical embedding 
separating a Minkowski from a black hole phase in the D4/D8 system, i.e. in the r-A plane. 
We will analyse this system separately in section |j. 

Let us now turn to the vector fluctuations on the probe branes. The equations gov- 
erning them follow from an expansion of the induced action in powers of the abelian field 
strength. At lowest order one has 



SVcctor OC / d xdAd0 6 _ p e 



9F MN F PQ g MP g N Q , 



(2.4) 



(and similar for the p-9 and r-A system) where g denotes the induced metric. The spatial 
components of the fields which describe the massive vector mesons arise from the expansion 



F 0a = VV) ( A ) > F ^ = H ^ W(n) W ■ 



(2.5) 



We will focus solely on thermal pole masses of the vector mesons, which are denned by 
-mtsBa ■ The equation of motion for the modes ip( n ) which then follows reads 



Ox 



i 2 -A 

+ m (n) e v 



(2.6) 
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Figure 2: The Schrodinger potential for Minkowski embeddings in the D3/D7 system, 
close to criticality, for various values r$ of the intersection with the vertical axis. The 
green dashed line indicates the asymptotic potential in the box, given by ( [3.17 ); the 
vertical dashed lines mark the edges of the plateau (which will be defined in section 3.3). 



These equations can be solved in a series expansion near the tip of the brane, and then 
extended numerically using standard shooting methods. In the Minkowski phase, the 
discrete mass spectrum arises after imposing that the fluctuations are regular in the IR 
and normalisable in the UV. In the black hole phase, the discrete quasinormal spectrum 
comes from incoming boundary conditions in the IR and normalisability in the UV. 

In order to get a better insight into the structure of the spectrum, it is useful to 
perform a coordinate transformation and a wave function rescaling which brings Q2.6| ) in 
Schrodinger form. Following the notation of [14], we write the equation in the general form 



r(A) 



-d\1p(n) 



m^(n)=0. (2.7) 



.5?(A) 

(The functions T and £ of course implicitly depend on the probe brane embedding func- 
tion r(A)). By introducing a new coordinate and a rescaled wave function according to 

dA 1 T fr , . 

d^=S' An)=^(n), - = Ve' (2 ' 8) 
the problem is transformed into a Schrodinger type problem, 

-d?i> in) +ViP in) =E n iP {n) , with V = ^, E n = m\ n) . (2.9) 

We will see how to approximate the potential V analytically in the following sections. 
However, to get an idea of what we are aiming for, let us first give the results of a numerical 
analysis, valid for the embeddings in the r-A plane. 

The resulting potential for the Minkowski embedding is displayed in figure || where tq 
is the position at which the probe brane intersects the vertical axis of figure |l[ The a 
coordinate is bounded on both sides. As the critical embedding is approached, the potential 
becomes more and more rectangular, while the box size (i.e. the range of a) grows, until it 
eventually becomes unbounded. 

The result for the black hole embedding is shown in figure ||, where xo = cos $o is 
related to the value of the angular coordinate at which the brane intersects the horizon, 
as shown in figure |l|. Now the a coordinate is bounded only in the UV direction, and the 
potential develops a clear step-shape. As the critical embedding is approached, the step 
becomes infinitely long. 
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Figure 3: The Schrodinger potential for black hole embeddings in the D3/D7 system, 
close to criticality, for various values of the angle xo of intersection with the horizon. 
Green and black dashed lines have the same meaning as in 0. 



These figures suggest that it should be possible to find a simple description of the 
potential very close to the critical embedding. We will do that in the next section. A 
somewhat similar story holds true for the D4/D8 embeddings, which we will analyse in 
more detail in section ||. 



3. The D3/D7 system 

3.1 Zero-temperature potentials 

Before giving our analytic results at finite temperature, it is instructive to recall the struc- 
ture of the Schrodinger potential at zero temperature. Although we mainly focus on the 
D3/D7 system (p = 3), we have kept part of the discussion general so the D4/D6 system 
can be obtained (by setting p = 4). At zero temperature, there are two types of embed- 
dings. One is the conformal embedding of the D7-brane (when the constituent quark mass 
m q is zero), and the other one is the non-conformal embedding, when m q 7^ 0. In the for- 
mer case the induced metric on the worldvolume of the brane is exactly that of AdSs x S 3 , 
while in the latter case the induced metric approaches this only in the UV, i.e. near the 
AdS boundary. The embedding is trivial in these cases, given by r = m q , which leads to 
the vector fluctuation equation [15] 4 

2 

la A (A 3 a A vv)) + (A2 + ( ^ Vv) = o. (s.i) 

The Schrodinger potentials for the m q = and m q > are easily computed to be, 
3 

V rn =0 (a) = —^, A = —a~ , 

4cr z 



V mq7 c (a) = m q 



1 3 
+ 



( tan 2 (m g o") + 



(3.2) 



A = m q t&n(m q <j) 



2 ' 4 V v.- ■-</-/ 1 tan 2 (m 9 cr)< 

In the conformal case the potential goes to a constant (zero) in the IR, while it diverges 
in the UV. In the non-conformal case, the potential diverges both in the UV and IR. The 



4 The curvature radius L has been set to one, hence the m q variable has dimension length. The coordi- 



nates used are r = ucosO and A = u sin 9 as in [15], but different from the coordinates of appendix A.l 
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box shape of the potential V mq ^o obviously leads to a discrete spectrum, in contrast to the 
continuum spectrum of the conformal case. If one replaces V mq7 Lo with an infinitely high 
box of size 7r/(2m g ), the spectrum asymptotically becomes m? \ = 4m^ n 2 , in agreement 
with the large-ra limit of the exact spectrum found in [15]. Hence, we see that while in the 
UV both branes are the same (i.e. AdS), with a diverging potential, the main difference 
between the two comes from the IR, i.e. from near the Poincare horizon. 

Let us now turn to the T / case. In this case the analytic form of the embed- 
dings of the branes is unknown, and so one would perhaps expect the analysis to be more 
complicated than above. However, the key element is that, for branes close to the critical 
embedding, only the near horizon region of the background is relevant, a region which is 
quite different from the T = background discussed above. Brane embeddings near the 
horizon were analysed by [16, 17]. As we will show below, when we get closer and closer 
to the critical embedding, the region in Schrodinger coordinate which corresponds to the 
near-horizon region becomes larger and larger, growing without bound. In the strict limit, 
we end up with an extremely simple, flat potential of infinite length, only corrected by 
finite-size effects whose relative effect scales to zero. The simple potential then allows us 
to make an exact statement about the spectrum. 

3.2 Brane embeddings in the Rindler approximation 

Having illustrated the behaviour of the Schrodinger potential with some numerical results 
and having sketched our approach, let us now provide an analytic argument which will lead 
to an exact potential in the critical limit. We first need to summarise the brane embeddings 
close to the horizon. As in [17], we write the metric of the 8 — p sphere as 

dQj_ p = d9 2 + sin 2 6dn 2 6 _ p + cos 2 Odip 2 . (3.3) 

We will consider the probe embedding extended in three x-directions, wrapping the 6 — p- 
sphere at fixed (p and having non-trivial 0(r). We now want to zoom in a region near the 
horizon, for a brane which is near its critical embedding. This means a region close to 
9 = 0, r = 7-q. Let us choose a parametrisation 



u = u T + nTz 2 , ^{(^j* > £ ={f)^ X - (3 ' 4) 
Then, considering the limit where 

u — ut ttTz y ( L \ 4 



« i , e = f — « i , (3.5) 

ut ut L \ut J 

we find Rindler space 5 

ds 2 = -(2^T) Vdi 2 + dz 2 + dy 2 + y 2 dQl_ p + dx\ + . . . , (3.6) 



5 The relation with the p, x variables of figure [l] is, for the D3/D7 system, given by p — 1 = z/L, 
X = 1 — ^y 2 /L 2 , valid near p w 1, x ~ lj see appendix |A.l| for details. 
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while the dilaton approaches a constant. Notice the horizon is now at z = 0. The Rindler 
approximation ceases to be valid at distances set by 



ut 



L p = 3; 



-^Rindler = \ — ™ = { , u (3-7) 

The equation describing the probe embedding y{z) in this region is 

zyy+{yy-(6-p)z)(l + y 2 ) =0. (3.8) 



An important aspect of ( |3.8j ) is its scaling symmetry, z — ► Az, y — > Ay. This means that 
changing z or y does not really influence the shape of the embedding but just rescales it. 
The obvious solution of this equation gives the critical embedding, 



Vcrit = \/Q-pz. (3.9) 

If the brane does not touch the horizon, there is a family of smooth Minkowski embeddings 
parametrised by zq, the distance of the tip of the brane to the horizon. The solution near 
z = zq reads 



2/Mink — Zq 



z \ \ Z Q 




(3.10) 



For z > zo one can find the asymptotic solution [16] which shows that the embedding 
approaches the critical one y cr it up to corrections which scale as a negative power of z. 
Taking the limit zq — > corresponds to approaching the critical embedding. 6 If the brane 
intersects the horizon, there is a family of black hole embeddings parametrised by yo, which 
near z = reads 



yBH = yo 



6 — p z 2 _ / / z 

l + —r L ^2+°\ — 



(3.11) 



Again, the parameter yo measures the deviation from the critical embedding, and taking 
the limit yo ^ corresponds to the limit in which the critical embedding is reached. In 
the other regime z S> zq the embedding again approaches the critical one [16]. 

3.3 Structure of the Schrodinger potentials 

Let us now analyse the Schrodinger potentials for the spectrum of vector meson fluctuations. 
The dominant part of the potential will turn out to be located in the Rindler region, so we 
will first focus on the embeddings described in the previous subsection. 7 We will focus on 
a simple subset of the possible excitations, namely those described (in the ^o = gauge) 
by Ai = e~ tujt ip(z). In Rindler coordinates, the equation for the fluctuation is 

d z (y 6 ~ p z(l + y 2 )-^) + ^(1 + y 2 )^V = , (3.12) 



6 However, note that the limit is subtle. Since the expansion parameter is (z — zo)/zo one cannot first 
take 20 — > 0, but rather one has to simultaneously scale 2 — > 0. This is similar to the situation we had for 
the T = embedding and the limit m q — > 0. 

7 Several other properties of holographic QCD are also fully determined by the Rindler region; see e.g. [18] . 
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Figure 4: The approximations used for the potentials close to the critical embedding, to 
be compared with the exact numerical results in figure | and §. 

where u = • Using (|3.8|) this equation can be written as 



2nT 

d 2 z 4> + + ^r^V = • (3.13) 

z z^ 

The IR boundary conditions to be imposed on the wave function tp are regularity in the 
Minkowski case and incoming boundary conditions in the black hole case, 

V>Mink = 1 + 0(z - zq) , i/i bk = z- ia (l + 0(z)). (3.14) 

It is useful to recast equation ( |3.12| ) in Schrddinger form, as we did for the full problem in 
the previous section. In terms of the notation used in [14] we now have 

r = ^(l + y 2 )i, Z = ^±fl, ~=y^, ^ = ^. (3.15) 

The new radial variable a is defined through da/dz = S, as before. With these redefini- 
tions, equation ( |3.12| ), describing fluctuations in the IR, can be rewritten as 

, 2 r xl/ 7 ~2 7 rr d 2 a E (6-p) z 2 6-p+((6-p)/2-l)(a 2 y) 2 

-d a ip+V m ip = uj ip , with Km - 



S 2y 2 l + (c5,y) 2 

(3.16) 

Using the explicit solutions for brane embeddings, described in the previous subsection, we 
can compute the potentials. They are plotted in figure for the p = 3 case. 



Let us first discuss the potential for the critical embedding (3.9). It turns out to be 
constant everywhere in the Rindler region, with a value given by 

v m _ ( 6 ~ p) 2 (o 17 \ 



The a coordinate is related to the z coordinate by a = ^J7 — p\og z, so the point z = is 
mapped to a = — oo while the edge of the Rindler region in the UV is mapped to a = +oo. 
The potential for the critical embedding is thus a constant on a = (—oo, +oo). 8 



8 If one repeats the computation for a scalar excitation y — y c i(z) + 5y(t,z), one gets for the critical 
embedding Krit = which 

is negative for the cases of interest p = 3, 4, signalling unambiguously 
the appearance of scalar tachyons in this channel near the critical embedding. These tachyonic modes were 
first explicitly discussed in [19]. 



- 9 - 



*" A 

Figure 5: Schematic indication of the three regions defined in the text. The light shaded 
area is the Rindler region. For a fixed embedding, the segment before entering the dark 
shaded area is the "deep IR" , the segment after exiting is the "UV" while the segment 
inside the dark shaded area is the "intermediate IR" . The latter grows, in the Schrodinger 
coordinate, as — 21og(r — l/\/2) or — log(l — Xo)- 

We have already seen numerical evidence for the formation of this plateau both from 
the Minkowski side as well as from the black hole side (figure || and ||). In the Minkowski 
case, the plateau arises from an approximately square well, while in the black hole case 
it arises from a step-like potential (note that the step is a generic feature of black hole 
embeddings [20]). We thus want to approximate these potentials as in figure [|. 9 Using the 
results of the present section, we can now make explicit how the lengths of the long flat 
segments of these potentials depend on the distance to the critical embedding. 

For all embeddings, three regions can be distinguished. First, there is a "deep IR 
region", where the near-critical Minkowski and near-critical black hole embeddings still 
deviate from the critical embedding. Somewhat further out, but still in the region where 
the Rindler approximation makes sense, there is an "intermediate IR region" or "plateau" . 
In this region the embeddings are well-approximated by the critical embedding. Finally, 
for z ~ -RRindien we leave the regime of validity of the Rindler approximation, and we have 
no analytic control anymore over the shape of the embeddings in this "UV region" . The 
three regions are shown in figure ||. 

The potential in the intermediate IR region is essentially constant, because the brane 
is there almost the same as the critical one. The key observation is that, as we argue below, 
the size of this intermediate region in the Schrodinger coordinate diverges logarithmically 
as the critical embedding is approached. The other two regions ("deep IR" and "UV 
region"), on the other hand, have a size which is independent of zq or i/q. Hence, as 
observed numerically, when we approach the critical embedding (either from the black hole 
or from the Minkowski side) the Schrodinger potential becomes constant, due to the fact 
that the size of the intermediate IR region in the Schrodinger coordinate diverges. 

9 For Minkowski embeddings in the D4/D6 case, the Schrodinger potential does not diverge at z = Zq 
but goes to a finite constant. However, regularity of the fluctuation requires ip to vanish at the finite value 
am — a(z — Zo)- Thus, a box potential is also a good approximation in this case. 
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Figure 6: Schrodinger potentials for the D3/D7 vector excitation in Rindler space, as 
well as an estimate of the left-hand side of the plateau, for the Minkowski and black hole 
embeddings respectively. 



Let us now compute the lengths of the various regions. First, in order to show that 
the deep IR region has a size which is independent of zq or yo> we make use of the scaling 
symmetry of ( |3.8| ), Suppose that we have a Minkowski solution 

y = fi(z) , which satisfies y{zo) = 0, (3.18) 

i.e. zq is the point where the brane is closest to the horizon. Using scaling symmetry, there 
is then another solution 

U = hiz) = /i ( z— ) , which satisfies y(zo) = 0. (3.19) 
V z oJ z 

Next, let us introduce a threshold A, which is a fixed and small positive number. When 
the distance between two branes is smaller than A, we will consider them to be coincident. 
Assume now that we have found a point z = zi e f t for which the distance between the first 
embedding and the critical one is equal to this A, 

h feft) " \/6 3 P zjdt = A • (3.20) 

This tells us immediately that, at a point z\ e f t = (z\ c { t /zQ) zq, the distance of the second 
embedding to the critical one satisfies 

ft \ ft \ t ( ^ left z 17, ?Left z A . A -r . - 

h{Z\ c { t ) - /crit(zieft) = fl Zq— — - ^6 - p Z ^- = — A < A It Z < Z . 

V Zq ZqJ Zq Zq Zq 

(3.21) 

In other words: if we define the end of the deep IR region to scale like linearly in zq, and if 
we fix the proportionality constant for some value of zq, then the approximation is better 
for all embeddings which are closer to critical. We can thus define the deep IR region to 
be zq < z < Azq for some A which is independent of zq. 

The length of the deep IR region in the Schrodinger coordinate a now comes out to 
be independent of zo, because d z y is a function of z/zq only (this is a consequence of the 
scaling symmetry). Therefore, using the value of E in ( p. 15 ), we find 



(1 + Y 2 (z/zq)) 1/2 ^_ (l + Y\q)) l/2 

z 



Ml 



A^ n p k IR =/ V ~ ' ' 7"""" dz= I v^— ^ d q, (3.22) 




Figure 7: Numerical values (for the D3/D7 system) for the width of the plateau in the 
Schrodinger a coordinate, for the Minkowski and black hole embeddings, together with 
their analytic approximations (3.23) and (3.24). 



where Y(z/zq) = d z y. This expression is independent of zq. The length of the "intermediate 
IR" region, Azo < z < ei?Ri n dier for some e < 1, can be computed using the critical 
embedding. 

We now concentrate on the D3/D7 case, for which the edges of the plateau are given 
by A ~ 3 and e ~ 1/5. One finds 

Aa^ e k au « f L ~dz = 2 log ^- « -21og(r - 1/V2) - 6.11 . (3.23) 

The numerical factors e and A are of course somewhat arbitrary, but the scaling with ro 
is fixed by the critical embedding. This scaling is reproduced well by the numerical data, 
see figure |7|a (the black dashed vertical lines in figure |2| indicate the edges of the plateau) . 
Finally, there is a UV region which is out of reach of the Rindler approximation. However, 
this UV region is similar to the one for zero-temperature embeddings, i.e. has finite size in 
the Schrodinger coordinate, independent of ro- 

A similar logic applies to the black hole embeddings. The agreement with the critical 
embedding (or equivalently, the start of the plateau) again occurs at around z 3yo, and 
the Rindler approximation breaks down at z ~ L/5. We then find that the plateau region 
scales as 

AogSeau ~ f ~ dz = 21 °S IT- ~ - " *0) " 6 - n ' ( 3 - 24 ) 

Again, this scaling with xo agrees well with the size of the plateau as measured from the 
numerical data, see figure ^b (the black dashed vertical lines in figure ^ indicate the edges 
of the plateau). There is also again a UV region, which looks like the UV part of the 
potential at zero-temperature, which has finite size in the Schrodinger coordinate. 

We thus see that in both cases, we can determine the size of the plateau as a function of 
the distance to the critical embedding. The plateau becomes of infinite length, and all other 
IR and UV features of the potentials get infinitely pushed away. They remain of finite size, 
however, and their contributions to the spectrum scale to zero in the critical limit (some 
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caricature potentials in which this can be seen explicitly are discussed in appendix A. 2 ). 
Therefore, we conclude that the mass spectrum can be approximated by using a square 
well potential on the Minkowski side, and a square step when one approaches from the 
black hole side (see figure ||) , where the plateau sizes scale as in ( |3.23| ) and ( |3,24 ) . 



3.4 Meson and quasi-normal modes 

Having obtained simple potentials for the two types of brane embeddings, let us now look at 
the spectrum. For the Minkowski embedding, the masses are simply given by the solution 
of the Schrodinger equation in a box of height V cr n ■ This yields 

y/Cbl - F cri t ~ — , n EN . (3.25) 

CUV — CIR 

However, the box size o"uv — °ir depends on zq, and diverges in the critical limit. Thus, all 
meson masses tend to \Acrit > as alluded to in the introduction. Verifying this numerically, 
on the other hand, is slightly subtle because the masses approach the limit only slowly. 

For the black hole embedding, we need to consider a potential with two regions. 
In the first region, with range — oo < a < cr s tep, the potential vanishes. The solution 
of the Schrodinger equation which has incoming boundary conditions at a = — oo is 
iprj\ = e~ luja . To simplify notation, we will set ouv = In the second region, with 



range c s tep < o" < 0, the potential has the value given in (3.17). The solution now reads 
■0(2) = Asin(\/Cd' 2 — Vcrit f) where A is a constant to be determined and we have imposed 
ip(a = 0) = since there is an infinite wall. By matching the values of ip and its derivative 
at a = (T s tep, we find a single equation for the values of the quasi- normal frequencies, 

iQ Sin(v / U> 2 - V crit CTstep) + V & 2 ~ Vcrit COs( \J OJ 2 - V CTit (T s tep) = . (3.26) 

Given Vcrit and <T s t ep , one can find numerically the complex values of u that solve this 
equation. However, let us look at a particularly interesting limit, 



k S tep|V^rit>l- (3-27) 



Since in our problem Vcrit is fixed to be the constant in ( [3.17 ), this is satisfied for large 



(negative) cr s tep- Remembering that |<7 s tep| grows unbounded when the flavour brane ap- 
proaches the critical embedding (jjq — ► 0), this is in fact the limit in which we are interested. 
But now we can obtain Co as an expansion in (o" s t ep VK;rit) _1 ) 



n 2 7r 2 . „ n 2 7r 2 



- VVcrit I ±1 ± - 2 - + iF + O (a step V^)" 4 ) | , n £ N . 

2a step V cnt V 

(3.28) 

Here F is the first coefficient in this expansion which depends on the finite-size details of 



the potential; for the step potential of figure f| we have F = 1 (for details see appendix A. 2 ). 

This expression is rather appealing. The real part is the one obtained if one had an 
IR infinite wall at <7 s t e p- That is, it corresponds to the meson masses of the Minkowski 
embedding if one identifies <T ste p with <tir there. In fact, from figures | and | we see that 
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Figure 8: Behaviour of the first quasi- normal mode of the D3/D7 system with a black 
hole embedding, as a function of the distance xo to the critical embedding. The blue 
dotted line represents the analytic result obtained by using a step potential (where F = 
0.84 was used in ( [3. 28 )). The values xo = 0.94 and xo = 0.9621 correspond to the points 
where the embedding becomes metastable and unstable, respectively. 



these values are comparable if zq and yo are of the same order. On top of that, there is a 
small, negative imaginary part (remember <7 s t ep < 0). 

The numerical analysis of the quasi-normal frequency spectrum away from the critical 
embedding, as presented in [4] for a scalar excitation, requires some care if one wants to 
extend it to the regime of near-critical embeddings. In the coordinates used here, we found 
reliable numerical results for the first quasi-normal mode, which match smoothly onto our 
analytic results. The result for the D3/D7 system is shown in figure ||. Equation ( 3.2S| ) 
thus provides us with an endpoint for the vector analog of the spiralling trajectories of the 
quasi-normal modes shown in figure 4 of [4]. 

Summarising the present section, we have seen that any finite number of meson masses 
in the D3/D7 system come down to a single value VKirit- Similarly, any finite number of 
quasi-normal modes come down to the same value as the critical embedding is approached 
from the black hole phase. Since the masses and quasi-normal modes are labelled by a 
single excitation number n, there is in principle an infinite number of ways to relate the 
states in these two spectra, and the connection between mesons and quasi-normal modes 
thus seems more subtle than previously suggested [5, 6]. Note that the degeneracy is not in 
contradiction with the Von Neumann & Wigner theorem [13]. In our case, a finite change 
of the parameter (here given by T/m q ) leads to an infinite change of the potential (the 
logarithmic growth which diverges at the critical embedding), while the Von Neumann & 
Wigner theorem requires the potential to change smoothly. 



4. The D4/D8 system 
4.1 Schrodinger potentials 

The D4/D8 embeddings used in the Sakai-Sugimoto model are different from the cases 
analysed so far. Instead of being described by a curve in the r — A plane, these embeddings 
are given by curve in the u— r subspace, i.e. these are more like the "equatorial embeddings" 
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of the D3/D7 system. There is a family of Minkowski- type embeddings, labelled by the 
lowest point uq of the D8 brane in the n-direction. The black hole embeddings are labelled 
by the distance in the r direction between the two halves of the D8 brane, but this distance 
is a modulus and the spectrum is independent of it. 

In the intermediate-temperature phase, where chiral symmetry is broken, the vector 
meson fluctuation equation reads [10] (we set ut = L = 1 from now on) 

du^T^Uuf'd^) +u-^f 4 (u)-^mf n) ^ n) = 0. (4.1) 

where 

7 = - . (4.2) 

u 8 h(u) - ug/4(«o) 

Following the by now familiar procedure we thus construct the Schrodinger radial variable a 
through 

a(u) = ± / ■ . (4.3) 

Juo V^lVu 5 (u 3 - 1) - u*(v% - 1) 

We have defined a = as the tip of the brane (u = uq) and the two signs correspond to the 
two branches of the brane worldvolume. For any uq > 1, it is clear that <7 max = a(u = oo) 
is finite, so a lives in a finite range. However, as uq — > 1, cr max diverges. The rate of 
divergence can be computed by making a coordinate change 5 = (n/no) 3 — 1, which leads 
to 

1 r°° (i + S)~h o 

^max(e) = , / ' ' d^ , where e = 1 - u . (4.4) 

" Jo Vovo + e 

This integral can be evaluated by splitting it into two regions, < 5 < \fe and yfe < 5 < oo. 
This leads to two integrals which can be evaluated analytically, up to corrections which 
vanish as y/t. The result is 

<w(e) = g (v^vr + 8 log 2 + log 3 - 2 \og{u - 1)) + 0{\fe) . (4.5) 

The Schrodinger potential in the u-variable is: 

_ n 5 (2u 6 + 2u 3 - 4) + u 5 (u 3 - l)(6n 3 - 9) 

Its limiting behaviour at the tip and far away from it is given by 

V(u = n ) = 8U ° ~ 13 5 U ° + 5 , y(n^oo) = oo. (4.7) 

The dependence of V on the a variable can be computed numerically. The result for 
various values of no is plotted in figure ^. These plots suggest that it should be possible to 
find two square well potentials, which provide an upper resp. lower bound on the potential. 



The lower bound is obtained by using a width which is twice the expression (4.5) (the red 
dashed outer box in figure |9|). 

For the width of the upper bound potential, we pick a point G{ on the a axis for 
which we can show that the potential goes to zero in the limit Uq — > 1, and which is such 
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-4-2 2 4 -4 -2 2 4 -4 -2 2 4 

Figure 9: The Schrodinger potential for the Sakai-Sugimoto system in the phase with 
broken chiral symmetry, for various decreasing values of the position uq of the tip of the 
brane. The dashed red lines (outer box) denote the potential used for the computation 
of the lower bound, while the dotted black lines (inner box) denote the potential for the 
upper bound (with q = 2). 



that cij — > oo in this limit. A suitable point is defined by 5 = e 1 /* 3 for some positive q. This 
point satisfies 



1 

3./„ 



■Mi 



1 



dS 



log 2 



Q 



3q 



1 g-i 
- loge + 0(e i 



(4.8) 



The potential V{(Ji) behaves as 0((uq — l) 1 ^ 9 ) and thus goes to zero. For sufficiently 
small no the potential is monotonic, and thus we conclude that V(a < <Ji) goes to zero as 
well. Furthermore, we can take the formal limit q — > oo in (|4.8| ). In this limit, the box 
size again diverges logarithmically as no —* 1, with the same coefficient of the logarithm as 
in ([4.5|). We will use a box with width twice that of (4.8) for q —* oo to obtain an upper 
bound to the potential (the black dashed inner box in figure ||) . 

In the phase where chiral symmetry is restored, the vector meson fluctuations are 
governed by the equation 



-u l l 2 {l - u-Z)d u (u b ' 2 {l - u- 3 )d u tf;{u)^ = u 2 i){u) . 



A Schrodinger equation is obtained by using 
1 1 



u 



V2(l 



u 



-3-] 



u; 



3 / 2 (l 



1/2 



ip = . 



(4.9) 



(4.10) 



The Schrodinger coordinate, obtained using da/du = S, is a complicated function of u 
which does not admit a simple inverse. We therefore restrict to a numerical plot of the 
potential in figure [TO. In the u coordinate it is given by 



V 



u° + u 6 



2w> 



(4.11) 



This potential is quite featureless and we will hence resort to a numerical determination of 
the quasi- normal mode spectrum for the D4/D8 case (this potential is similar to the one 
for the equatorial D3/D7 embedding). 
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Figure 10: The Schrodinger potential for the Sakai-Sugimoto system in the phase where 
chiral symmetry is restored. The a coordinate spans the half-line — oo < a < 0. 



4.2 Mesons and quasi-normal modes 

Using the square well potentials which we have shown in the previous section to be upper 
resp. lower bounds for the full Sakai-Sugimoto potential, it is easy to compute the mass 
spectrum in the uq — > 1 limit. From the box with a width given by twice (4.5) we find 



nir 



m {n) {u ) > —— v , nGN, (4.12) 

2 ( + | log 2 + I log3 - 1 log(n - 1)' 



whereas the box with width given by twice (4.8) we get 



nir 



2 (I log 2 - I log 3 - I log(u - 1)) 
To re-instate the L and ut dependence, note that 

m = , m = . (4-14) 



u 



1/2 A7tT 



T 



These expressions both go to zero as uq — > 1; see figure |ll] for a comparison with the 
numerically determined spectrum. 

Because the potential in the high-temperature phase ( 4.1 1| ) does not take a simple form, 
we again have to resort to numerics in order to find the quasi-normal mode spectrum. The 
first six vector meson masses are given in figure (^) (each of these is twofold degenerate 
because of chiral symmetry). A fit to a linear function yields 

Q n = ±(0.31 + 1.305n) - (0.145 + 0.752n)i , (4.15) 

(where u is related to to by a relation similar to ( [4. 14| ) ) . The numbers agree with the 
analysis presented in [21]. 

Comparing the real parts of the quasi-normal modes with the meson masses in the 
intermediate-temperature regime, we see that there is no room to connect the two in a 
smooth way. The masses of any finite number of mesons come down to zero as the critical, 
overheated embedding is approached. In contrast to this, the quasi-normal modes take 
values which are equally spaced in the complex plane, with a lower bound on the real part 
which is non-zero. An identification along the lines of [5, 6] thus does not seem possible 
for the Sakai-Sugimoto model. 
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n 






1 


±1.61 


- 0.90i 


2 


±2.93 


- 1.65z 


3 


±4.23 


- 2.m 


4 


±5.54 


- 3.16? 


5 


±6.84 


- 3.91s 


6 


±8.13 


- 4.66i 



Figure 12: The spectrum of the first twelve quasi-normal vector modes of the chirally 
symmetric phase of the Sakai-Sugimoto model (each mode is twofold degenerate because 
of chiral symmetry) . 

5. Conclusions and discussion 

We have analysed the form of the Schrodinger potential for vector meson fluctuations in 
various holographic duals to gauge theories with matter. These potentials exhibit some 
generic structure which we expect to see for more general classes of models. In the Rindler 
region, a Minkowski embedding will lead to a box-shaped potential, while a black hole 
embedding will lead to a step. When no Rindler approximation can be made, which is the 
case for the equatorial embeddings, more of the background can be seen, and the potential 
is typically more complicated. 

In the non-equatorial type transition of the D3/D7 system, we have shown how any fi- 
nite number of meson masses come down to a fixed value. Similarly, when the embedding is 
approached from the black hole side, any finite number of quasi-normal modes approaches 
this same value. Lacking any quantum numbers to distinguish them, the connection be- 
tween these spectra is thus more subtle than previously suggested. 10 



10 In [5] it was observed that peaks appear in the spectral function, which is consistent with our results. 
However, our calculations show that these peaks all converge onii = 0.75. The spectra are symmetric only 
for the artificial choice of fixed zo/yo- 
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For the D4/D8 embeddings ("equatorial"), we have shown that all meson masses come 
down to zero. This is consistent with the fact that there is a proper chiral symmetry 
restoration after the transition, in sharp contrast to the D3/D7 system. The quasi- normal 
mode spectrum, on the other hand, is equally-spaced with a non-zero lowest frequency 
mode. This result does not support the idea that quasi-normal modes are continuously 
connected to meson modes. It is conceivable that, if one would introduce a bare quark 
mass into the Sakai-Sugimoto model (along the lines of [22-24]), its behaviour would be 
more similar to the D3/D7 case. It would also be interesting to verify if the qualitative 
results we have obtained for the D4/D8 model are realised in its non-critical version [25], 
whose thermal phase structure was analysed in [26]. 

We have not yet touched the issue of chemical potentials. In the presence of a baryon 
chemical potential, the spectral function exhibits, for large enough m q /T, sharp peaks close 
to the critical embedding [6]. Moreover, these peaks converge on the meson masses for the 
Minkowski embeddings. However, while the embeddings for large chemical potential are 
very close to the Minkowski embeddings (and hence the agreement in the spectrum), they 
are not related to each other by a small deformation in the gauge theory parameter space. 
The situation is reminiscent of a strong/weak duality in chemical potential. It would be 
interesting to understand this behaviour for isospin potentials as well [27]. 
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A. Appendix 

A.l Coordinate frames 

We collect here a number of technical details related to the coordinate choices made in the 
main text. We will rewrite the metric (|2.l| ) in two different ways, as depicted in figure |l|. 
One coordinate frame is better adapted for computations with Minkowski embeddings, 
while the other one is more suitable for black hole embeddings. 

The first step is to define coordinates such that the second part of the metric ( |2.1| ) 
becomes conformally flat. This is achieved with 



2 



P-3 

u 2 




(A.l) 



The metric fl2.1| ) becomes 




f p dt 2 + 5ijdx 



) 




(A.2) 



We can reexpress the flat factor inside the last bracket by defining p 2 = r 2 + A 2 , and we 

get 

7-p 

ds 2 = (I) 2 (- fpdt 2 + tfydxWj + L^KpfdX 2 + X 2 dn 2 _ p + dr 2 + r 2 dc/> 2 ) . (A.3) 

This is the r — A system depicted in figure p], an d for p = 4, it corresponds to the coordinate 
system defined in [2]. We have found this frame convenient to numerically analyse the 
Minkowski embeddings which reach A = at some r = ro, which parametrises the one- 
parameter family of embeddings (Schrodinger potentials for various values of ro in the 

2 

D3/D7 case are plotted in figure ||). The critical embedding corresponds to ro = 2p~ 7 . 

We now define a p — x coordinate system (as shown in figure |l|) following [17]. We first 
rescale the p coordinate defined above as: 

p = 2^p, (A.4) 

such that the horizon is at p = 1 and write the metric of the sphere as dOg_ p = d6 2 + 
sin 2 8dfl6-p + cos 2 9d<p 2 . If we now define x = cos we get 

7-p 

ds 2 = (^) 2 (- f P dt 2 + <%dxW) 

+ ^rf(^ + T ^ + (l- X 2 )d^ p + xW) , (A.5) 

where u should be understood to be a function of p. This coordinate system is convenient 
to study the black hole embeddings, parametrised by xo = x\p=l (with Xo = 1 being the 
critical embedding). Some Schrodinger potentials for various values of xo i n the D3/D7 case 
are shown in figure ||. In this set of coordinates, Minkowski embeddings are parametrised 
by Po — p\x=i- For embeddings near the critical one, there are simple relations between 



these parameters and t/o,zo as defined in section ^2. For the D3/D7 case (p = n = 3), 
they read: 



y = V2 y/1 - xo L (1 + 0(1 - xo)) z = L(p -l) (l + O( P0 -l)) (A.6) 

These relations, together with po = \/2ro, were used to obtain the last equality in equa- 
tions JSP!! ) and dj| ). 



A. 2 Quasinormal frequencies in toy potentials 

We have argued in the main text that for the D3/D7 and D4/D6 cases near the crit- 
ical embedding, the Schrodinger potential presents a plateau of diverging width in the 
a-coordinate. At its edges, there are some IR and UV features which remain of finite size. 
Since the plateau is infinitely big as compared to these finite size details, one expects that 
these details only affect the quasinormal frequencies (or meson masses) in a subleading 
way. In this appendix, we will explicitly check this for two toy Schrodinger potentials that 
approximate the near critical Schrodinger potential for D3/D7 black hole embeddings. 
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By looking at the plots in figure [| one sees that there is a bump before the plateau and 
a dip after the plateau, before the infinite wall. We stress that the size of the bump and 
dip remain constant as the plateau grows unbounded when the embedding approaches the 
critical one. Let us, accordingly, consider the following approximation for the Schrodinger 
potential which we insert in (2.9): 



V = for a < dstcp + x 1 , 

V = V CTit + hi for CJ ste p + 2i < O < (T stcp , 

V = V CTit for a stC p < cr < x 2 , 

V = Vcrit - h 2 for x 2 < a < , 



(A.7) 



and an infinite wall at a = 0. We want to consider V^-it, h±,h 2 > and a s t ep , x±, x 2 < 0, 
in a limit <7 s t C p — > — oo with V CT n,hi,h 2 ,xi,x 2 fixed. We can obtain a solution of the 
Schrodinger problem with IR incoming boundary conditions (if) = e~ luJa in region 1) as a 
series in l , . The solution is the same as in the simple step potential, see (|3,2§| ), 



^"stcpV V^ r it 

with F depending on the various parameters 



F 



hih 2 V CI it cosh(\/ h\x\) — ih\ \fh% sinh(y / h\X\ i 
cosh(v / /ii2i) \/ hih 2 V CT - lt x 2 + i\/hih 2 V CT i t - y/h\V cr i t ta^y^a^ 



(Ai 



+ sinh(A//iixl) -ihi\f h 2 V cr - lt X2 - \fh 2 V cr it + ih\ \A^rit tan(y / /^x 2 



If we insert h\ = h 2 = .1, xi = x 2 = —2, which is approximately appropriate for the 
potentials we are dealing with, and also use V cr \t = 9/16, we get: F « .65 — 1.72i. 

One may also worry if the exponential tail of the Schrodinger potentials, which is 
proportional to e 2fT as a — ► — oo could affect the computation. Let us address this point by 
considering another toy potential slightly more complicated than the one in figure |j: 



V = V ciit e 

V = V cr[t 



2(7- 2(7; 



step 



for a < CTstep , 

for CJstep < cr < , 



(A.9) 



Again, one can make an expansion for the o" s tep - * —oo limit and the solution for the 
quasinormal frequencies is again JOaP where F now can be written in terms of the modified 
Bessel function: 

-2i I , /Tr-{y/V~ t ) 



If we insert Vcrit = 9/16, we find F ~ 1.1 — .321 

These two examples explicitly show how in the cr s tep — ► — cxd limit, the details of 
the potential only appear in uji n \ at order (a s t e p\Acrit) -3 (which is however the leading 
imaginary part for the quasinormal frequencies) . In summary, this shows that finite details 
do not change the discussion based on the infinite plateau. 



- 21 - 



References 



[1] J. Babington, J. Erdmenger, N. J. Evans, Z. Guralnik, and I. Kirsch, "Chiral 
symmetry breaking and pions in non-supersymmetric gauge/gravity duals", Phys. 



Rev. D69 (2004) 066007, hep-th/0306018 



[2] M. Kruczenski, D. Mateos, R. C. Myers, and D. J. Winters, "Towards a holographic 
dual of large-A c QCD", JEEP 05 (2004) 041, |hep-th/031127C . 

[3] O. Aharony, J. Sonnenschein, and S. Yankielowicz, "A holographic model of 

deconfinement and chiral symmetry restoration", Ann. Phys. 322 (2007) 1420-1443, 
nep-th/0604161 . 

[4] C. Hoyos-Badajoz, K. Landsteiner, and S. Montero, "Holographic meson melting", 



JEEP 04 (2007) 031, hep-th/0612169 



[5] R. C. Myers, A. O. Starinets, and R. M. Thomson, "Holographic spectral functions 
and diffusion constants for fundamental matter", JEEP 11 (2007) 091, 
|arXiv: 0706 .0162 [hep-th] . 



[6] J. Erdmenger, M. Kaminski, and F. Rust, "Holographic vector mesons from spectral 



functions at finite baryon or isospin density" , |arXiv : 0710 . 0334 [hep-th] 



[7] K. Peeters, J. Sonnenschein, and M. Zamaklar, "Holographic decays of large-spin 



mesons" , JEEP 02 (2006) 009, |hep-th/05 11044 



[8] QCD-TARO Collaboration, P. de Forcrand et al, "Meson correlators in finite 



temperature lattice QCD", Phys. Rev. D63 (2001) 054501, liep-lat/0008005 



[9] M. B. Oktay, M. J. Peardon, J. I. Skullerud, G. Aarts, and C. R. Allton, 

"Charmonium properties in the quark-gluon plasma" , PoS LAT2007 (2007) 227, 
|arXiv: 0710 .2795 [nep-lat] . 



[10] K. Peeters, J. Sonnenschein, and M. Zamaklar, "Holographic melting and related 
properties of mesons in a quark gluon plasma", Phys. Rev. D74 (2006) 106008, 
lhep-th/0606195 . 

[11] F. Karsch, "Lattice QCD at finite temperature and density", Nucl. Phys. Proc. 
Suppl. 83 (2000) 14-23, |hep-lat/9909006| . 

[12] S. A. Gottlieb et al., "Thermodynamics of lattice QCD with two light quark flavours 
on a 16 3 x 8 lattice. II", Phys. Rev. D55 (1997) 6852-6860, |hep-lat/9612020 . 

[13] J. von Neumann and E. Wigner, "On the behavior of eigenvalues in adiabatic 
processes", Z. Phys. 30 (1929) 467-470. 

[14] E. Schreiber, "Excited mesons and quantization of string endpoints", 
lhep-th/0403226 . 



- 22 - 



[15] 
[16] 

[17] 
[18] 
[19] 
[20] 
[21] 
[22] 
[23] 
[24] 
[25] 

[26] 
[27] 



M. Kruczenski, D. Mateos, R. C. Myers, and D. J. Winters, "Meson spectroscopy in 



AdS/CFT with flavour", JEEP 07 (2003) 049, |hep-th/030403; 



V. P. Frolov, A. L. Larsen, and M. Christensen, "Domain wall interacting with a 
black hole: A new example of critical phenomena", Phys. Rev. D59 (1999) 125008, 
|hep-th/9811148| . 

D. Mateos, R. C. Myers, and R. M. Thomson, "Holographic phase transitions with 
fundamental matter", Phys. Rev. Lett. 97 (2006) 091601, |hep-th/0605046| . 



K. Peeters and M. Zamaklar, "Dissociation by acceleration", JEEP 01 (2007) 038, 
|arXiv: 071 1.3446 [hep-th] . 



D. Mateos, R. C. Myers, and R. M. Thomson, "Thermodynamics of the brane", 



JEEP 05 (2007) 067, |hep-th/0701132 



G. T. Horowitz and V. E. Hubeny, "Quasinormal modes of AdS black holes and the 
approach to thermal equilibrium" , Phys. Rev. D62 (2000) 024027, |hep-th/9909056| . 



N. Evans and E. Threlfall, "Mesonic quasinormal modes of the Sakai-Sugimoto 



model at high temperature" , arXiv : 0802 . 0775 [hep-th] 



O. Bergman, S. Seki, and J. Sonnenschein, "Quark mass and condensate in HQCD", 



JEEP 12 (2007) 037, |arXiv: 0708 .2839 [hep-th] 



A. Dhar and P. Nag, "Sakai-Sugimoto model, Tachyon Condensation and Chiral 
symmetry Breaking", JEEP 01 (2008) 055, |arXiv : 0708 . 3233 [hep-th] . 



R. Casero, E. Kiritsis, and A. Paredes, "Chiral symmetry breaking as open string 



tachyon condensation", Nucl. Phys. B787 (2007) 98-134, tiep-th/0702155 



R. Casero, A. Paredes, and J. Sonnenschein, "Fundamental matter, meson 
spectroscopy and non-critical string/gauge duality", JEEP 01 (2006) 127, 
|hep-th/051011C . 



V. Mazu and J. Sonnenschein, "Non critical holographic models of the thermal 
phases of QCD" , |arXiv: 0711 . 4273 [hep-th] . 



O. Aharony, K. Peeters, J. Sonnenschein, and M. Zamaklar, "Rho meson 
condensation at finite isospin chemical potential in a holographic model for QCD", 
JEEP 082 (2007) 1007, |arXiv: 0709 .3948 [hep-th]" , 



- 23 - 



